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Abstract 

Numerical renormalization group and conformal field theory work indicate 
that the two impurity Kondo Hamiltonian has a non-Fermi liquid critical 
point separating the Kondo-screening phase from the inter-impurity singlet 
phase when particle-hole (P-H) symmetry is maintained. We clarify the cir- 
cumstances under which this critical point occurs, pointing out that there are 
two types of P-H symmetry. Only one of them guarantees the occurance of 
the critical point. Much of the previous numerical work was done on models 
with the other type of P-H symmetry. We analyse this critical point using 
the boundary conformal field theory technique. The finite-size spectrum is 
presented in detail and compared with about 50 energy levels obtained us- 
ing the numerical renormalization group. Various Green's functions, general 
renormalization group behaviour, and a hidden SO {7) are analysed. 

75.20.Hr 
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I. INTRODUCTION 



While the single-impurity Kondo problem is by now rather well understood, the 
Kondo lattice model, of possible relevance to heavy fermion materials, presents additional 
difficulties.0 In particular there are two competing tendencies; the Kondo effect leads to 
a magnetic screening of the spins by conduction electrons while the RKKY interactions 
between the spins may lead to antiferromagnetism. The two-impurity Kondo problemiHB 
provides a simple model in which to study these competing effects. Of course, in this model, 
true antiferromagnetism with static moments (i.e. spontaneous breaking of rotational sym- 
metry) cannot occur. Rather the RKKY couplings tend to promote singlet formation by 
the two impurities. 

A very general form of the Hamiltonian has two s = 1/2 spins, symmetrically located 
about the origin, interacting with a Fermi gas. The Hamiltonian is: 

H-fiN = Ho + HK + H^^^^ (1.1) 

where Hq, Hk and -ffggjf are the kinetic energy, Kondo interaction and impurity self inter- 
actions respectively: 



^self = ■ S2 (1.2) 

v{k) is related to an Anderson model hybridization matrix element and K is the inter- 
impurity coupling. 

The physics of this model becomes simple in the limits K —>■ ±oo. When K — ^ +oo the 
two impurities lock into a singlet state. Therefore the Kondo coupling has no effect and the 
electron gas is completely unaffected by the impurities. Conversely, when K ^ —oo, the 
impurities become an effective single s = 1 impurity. As we shall see in the next section, 
the extended nature of this s = 1 impurity implies that two "channels" of conduction 
electrons interact with it. The two channel, s = 1 Kondo problem was studied by Nozieres 
and Blandin.! They concluded that, at low temperatures, the impurity gets magnetically 
screened by the conduction electrons. The remaining low energy conduction electron degrees 
of freedom are decoupled from the impurity but experience a ^ phase shift in both channels. 
This corresponds to a local Fermi liquid fixed point. The low energy electronic degrees of 
freedom are free electron- like; the many body interactions induced by the Kondo interaction 
lead to a simple phase shift at low energies. 

The behaviour of the system at intermediate values of K is less obvious. In general, 
one might expect that a local Fermi liquid description holds at low temperatures for all K. 
The phase shifts could vary continuously with K. However, it was argued by Millis, Kotliar 
and Jone^l that this cannot happen if P-H symmetry is maintained while K is varied. In 
this case, the zero-energy phase shift, if it is well-defined, can only be or 7r/2. Since the 
K — > ±oo fixed points are stable, it follows that there must be at least one point in the 
phase diagram not belonging to either phase, corresponding to some sort of phase transition. 
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Such a transition could, in principle, be first order. Alternatively, if it is continuous, the 
critical point cannot be characterizable by phase shifts. Numerical renormalization group 
(NRG) work indicated that the transition is indeed continuous. The critical point is of 
non-Fermi liquid type. On the other hand, quantum Monte Carlo (QMC) work has not seen 
such a critical point.0 NRG work on the Anderson model0 has seen the critical point only 
in the case of "energy-independent coupling constants". In both cases P-H symmetry was 
maintained. 

In the next section we re-examine the argument for a critical point in more detail. We 
show that there are actually two quite different types of P-H symmetry which may occur in 
models of this type; this is connected with the parity symmetry of the models. Only one of 
these P-H symmetries guarantees a phase transition. Much of the numerical work in Refs. 
( was done using models with the "wrong" type of P-H symmetry. 

Two of us have recently developed a new boundary conformal field theory technique 
to study such non-trivial critical points.0'iJii The essence of this method is that the local 
interactions can be replaced by a conformally invariant boundary condition^ in the low 
energy effective Hamiltonian. One of the main purposes of this paper is to explain the 
application of the technique to this problem and elucidate its various predictions. A brief 
presentation of some of these results was made earlier .0 In Sec. HI we review the reduction 
of the problem to a one-dimensional field theory and then discuss a convenient and some- 
what unconventional "bosonization" of the model in terms of Wess-Zumino-Witten matrix 
fields and Ising model fields. In Section IV we solve for the non-trivial critical point by 
hypothesizing the corresponding boundary conditions (there are only a few possibilities). 
We calculate the resulting finite size spectrum using this boundary condition and compare 
it to the NRG results, obtaining excellent agreement. In Sec. V we discuss the various 
scaling operators at the non-trivial critical point, the stability of the fixed point (including 
the effects of P-H symmetry breaking) and the behaviour of various Green's functions. The 
non-trivial fixed point possesses a remarkable hidden SO (7) symmetry.0 We explain this 
in detail in Sec. VI and comment on the possible connection with a physical picture of 
the 2-impurity Kondo problem based on abelian bosonization and "refermionization" .IllS 
In Sec. VII we summarize the disagreement with other calculations and the prospects for 
resolving them. 



II. EXISTENCE OF THE NON-TRIVIAL FIXED POINT AND PARTICLE-HOLE 

SYMMETRY 

We begin with a few comments on the Hamiltonian of Eq. ( p..2| ). Often, one does not 
consider an explicit inter-impurity interaction but only the indirect RKKY interaction, of 
second and higher order in J. More generally an additional interaction could be generated 
by other exchange processes not involving the conduction electrons. Often one considers a 
5-function Kondo interaction, with the impurities at ±R/2 in which case: 

v{k) = t;oe'^-^/2_ (2.1) 

The more general form that we consider allows for longer range hybridization. Note that, if 
e{—k) = e{k), the Hamiltonian is invariant under the parity transformation: 
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5"! <-> 5*2 

i^k ^ ^-h (2-2) 

corresponding to a reflection about a point in position space (the midpoint between the 2 
impurities.) For a tight-binding model, the origin of parity symmetry may not correspond 
to a lattice point. Nonetheless, we define our Fourier transforms with respect to it. 

For certain choices of the dispersion relation, e(fc) and hybridization matrix element, 
v(/c), the Hamiltonian of Eq. (1^) will have particle-hole (P-H) symmetry. We consider 
invariance of H under a general P-H transformation of the form: 

S\ — * 5*1 

5^2 ^ ^2. (2.3) 
where /c' is some function of k. Invariance of ifo requires: 

6(fc) = -e(fc'). (2.4) 

Invariance of the Kondo interaction requires (after a possible phase redefinition of v which 
doesn't effect Hk): 

v{k') = v{k)* 

v{-k') = v{-kye'" (2.5) 

where a is a ^-independent phase. For instance, for a cubic lattice tight-binding model at 
half-filling, 

k' = kQ- k, (2.6) 

where /cq is the nesting wave- vector: 

— * 

fco = (vr/a, 7r/a, 7r/a). (2-7) 



Hence, for the 5-function Kondo interaction of Eq. (|2.1| ), and appropriate choice of the 
phase of Vq] Eq. ( p.5|) is obeyed with: 

a = ko-R. (2.8) 

— * 

Assuming that R (the vector connecting the two impurities) is a lattice vector, the phase a 
can take on the values or vr. 

When particle-hole symmetry is present, there is also an exact SU{2) "isospin" symmetry, 
in addition to the normal spin symmetry.El [For a general discussion of this see Sec, II of 
Ref. ( |19D.] One of these isospin symmetry generators is simply the total electron number, 
corresponding to the 2;-component of isospin, P. The lowering operator, /~ is proportional 
to: 
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/ ^'%T^^'a> (2-9) 

— * — » 

where k' is the function of k occuring in the particle-hole transformation. 

The Hamiltonian of Eq. (|1.2|) can be reduced exactly to a one- dimensional one. To do 
this we define two fields: 

= I d^kS[e{k) - E]v{±k)i;^. (2.10) 

Note that, of the infinite set of fields of a given energy E, (corresponding to all points on 
the constant energy surface in /c-space) only these two appear in the Kondo interaction. We 
may define parity even and odd orthonormal linear combinations of these two fields: 

lPe,E = {ij+,E + i^-,E)/Ne{E) 

^o,E = {^+,E - tlj-,E)/No{E), (2.11) 

where: 

Ne,o{E) = J Sk5{E - e{k)]\v{k) ± v{-k)\^. (2.12) 

The anti-commutators are normalized as follows: 

{iPe.'iP%,} = 5{E-E'). (2.13) 

We may complete ifje and ipo into a complete orthonormal basis for each energy, E\ only 
these two fields appear in the Kondo interaction. Discarding the additional fields which 
decouple, we obtain an exact one-dimensional rewriting of the original Hamiltonian, with: 

Ho = j dEE[tPlj,^,^E + ^1,e^o,e] 

Hk = J dEdE'^^ [N,{E)N,{E')ijlEa^e,E' + iV„(E)iV„(E')<£;<?V^o,i?') ■ {Si + S2) 

+N,{E)N,{E') {iPIe^^cE' + 4,E^^e,E') ■ {Si - S2)}. (2.14) 

If the original problem had the particle- hole symmetry of Eq. (2^), then the one- 
dimensional problem also has P-H symmetry. The transformation of the fields ipe^E and ijjo,E 
can be deduced from: 

1p+,E 1pi,-E 

i,.^E ^ e*"^l,-i5, (2.15) 



which in turn follows from Eq. (|2.10| ) using Eq. p.3| )-( p75[ ). 



[We suppress the transformation of spin indices; it is the same as in Eq. (|2.3| ).] In the 
case a = we have [ using Eq. ( p^ ) and (p.5|)]: 

N,{-E) = N,{E) 

No{-E) = No{E) (2.16) 
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and 



i^o,E ^ ^l,-E- (2.17) 

N,{-E) = No{E), (2.18) 

i^e,E i^l-E 

^o,E^4,-E- (2.19) 

It can be directly verified that the P-H transformation of Eq. ( p.l7| ) or ( p.l9| ) are symmetries 
of the one-dimensional Hamiltonian of Eq. ( 2.14| ) if Ne{E) and No{E) satisfy Eq. ( p.l6| ) or 



For a = TT, we find: 



and 



(^TT^) respectively. 



It turns out that the model behaves very differently depending on which type of P-H 
symmetry is present. With the first type, a phase transition must separate the Kondo- 
screened and inter-impurity singlet phases; with the second type no transition is required 
or expected to occur. We reiterate the argument for a transition with the first type of 
P-H symmetry.! In a Fermi liquid phase, it should be possible to characterize the zero- 
temperature fixed point by phase shifts for the even and odd channels at -E = 0. These 
amount to boundary conditions relating incoming and outgoing operators: 

^out ^ ^2...^in 

^out ^ ^2..„^in _ (2.20) 

The Hermitean conjugate fields, iplo obey the same conditions with 5e,o —Se,o- Hence, 
if the first type of P-H symmetry is obeyed, both 6e and 6o can only take the values 
or 7r/2 [from Eq. ( p. 17 ) at E = 0, right at the Fermi energy]. [Note that the boundary 



conditions of Eq. ( p.20| ) only depend on 26 (mod 2tt) so that 5 = or vr are equivalent 
as are 6 = ±7r/2.] 5e = 5o = in the inter-impurity singlet phase and Se = So = 7^/2 in 
the Kondo-screening phase. Furthermore, both of these fixed points are absolutely stable 
(no relevant or marginal operators.) It then follows that there must be some sort of phase 
transition separating these two phases. It could be first order, or correspond to a non-Fermi 
liquid critical point at which the phase shifts are not defined. On the other hand, if we 
have the second type of P-H symmetry, 6e and 6o may take on arbitrary values subject 
only to the condition 6e = —So- Thus they may vary continuously from (0, 0) in the limit 
of cxo antiferromagnetic inter- impurity coupling to (7r/2, — 7r/2) in the opposite limit of oo 
ferromagnetic inter-impurity coupling. No transition is necessary in this case. 

The same argument can be made by considering possible potential scattering terms in 
the Hamiltonian. While the original Hamiltonian has no such terms, they will be generated 
in higher orders of perturbation theory. In general these take the form: 

Hps = J dEdE'[Ve{E, E')^{e^,^e' + Vo{E, E')^I,^^o,i?']. (2.21) 



6 



The phase shifts at zero energy contain terms proportional to ^£(0,0) and K(0, 0) respec- 
tively, for weak potential scattering, ignoring the Kondo interactions. With the first type of 
P-H symmetry: 

V.^oiE, E') = -V,^o{-E', ~E). (2.22) 

This implies K,o(0, 0) = 0; no phase shift at zero energy. On the other hand the second type 
of P-H symmetry implies: 

V,{E,E') = -Vo{-E',-E). (2.23) 

At zero energy this gives: Ve(0,0) = — Vo(0,0), allowing (equal and opposite) phase shifts 
at zero energy. 

In Sec. IV we explicitly study the stability of the non-trivial critical point under potential 
scattering. Our analysis shows that Ve(0,0) = — K(0,0) 7^ is a relevant perturbation. 



III. ONE DIMENSIONAL BOSE-ISING REPRESENTATION 

In this section we derive a representation of the two impurity Kondo problem which is 
convenient for studying the non-trivial critical point. The one-dimensional Hamiltonian of 



Eq. ( 2.14| ) is treated using a generalized bosonization method based on the Goddard-Kent- 
OliveB" coset construction. We end up representing the fermions by three Wess-Zumino- 
Witten non-linear a-models^ together with an Ising model. 



Taylor-expanding Ne,o{E) in Eq. (|2.14|) around E = 0, we obtain the non-derivative 
interactions: 

Hk = ^ J dEdE'l^ [jei)l^E^i)e,E' + Joi^l^E^i'o,E') " {Si + S2) 

+ Jm {i)l^E^i)o,E' + ^o,i5<T^e,E') " (5^1 " 5^2)| (3.1) 

where 

J, = 2N,{Q)\ Jo = 2No{Q)\ Jra = 2No{'d)No{Q). (3.2) 

Sums over spin indices are implied. We effectively obtain a two channel one-dimensional 
model. The origin of the two channels is simply the fact that the two impurities are at 
different spatial locations. In general, with n impurities at different locations, we would 
obtain an n-channel one-dimensional model at long wavelengths. 

The higher terms in the Taylor expansion of Ne,o{E) give various derivative interactions. 
These are irrelevant at the weak coupling (unstable) critical point, assuming the first type 
of particle-hole symmetry discussed in Sec. II. That is to say, they are irrelevant assuming 
they don't generate any potential scattering terms; this is the case with the first type of 
P-H symmetry. From our analysis of the non-trivial critical point we will conclude that 
these extra terms are also irrelevant at that critical point, assuming the first type of P-H 
symmetry. 
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We now see that, in the hmit of an infinite ferromagnetic inter-impurity interaction K, 
the model reduces to a single-impurity two-channel s = 1 Kondo problem, with different 
couplings, Je and Jo to the two channels. This difference is known to be irrelevant!. 

For analysing the critical point, we actually find it somewhat more convenient to intro- 
duce a different orthonormal basis: 



^Pl,2 = {ipe ± V'o)/V2 (3.3) 
The kinetic energy remains diagonal, as in Eq. ( |2.14| ) and the Kondo interaction becomes: 

Hk = -j dEdE'{J+[ipli^ai!i^E' + '4^1,e^'^2M ■ [Si + S2] + JnM^l^E^'^i^E' - 

^iE^^2,E'] ■ [Si - S2] + J^[i)\^E^i)2^E' + ^2,E^'4^lM ' ['S^i + ^2]} (3.4) 

where 

J± = ( Je ± Jo)/2 (3.5) 

Note that il)i,2 are not, in general, the same as il)± since A'^e 7^ No- We emphasize that this 
is the most general Hamiltonian consistent with particle-hole symmetry, up to operators 
which are irrelevant at the (unstable) zero-coupling fixed point, ie. in ordinary perturbation 
theory. 

The conserved isospin operators referred to above take simple forms in the effective 
one-dimensional theory: 



If 1 

= -j dE[lpl'^j^^o,l,E + ^'^l%i^a2,E] 
I' = j dE[i)i^i^E^ii_E + V'T2,i?V'i2 -b] (3.6) 

We note that when J_ = there is even more symmetry. Now the charges of the 1 and 
2 species of fermions are separately conserved and in fact we have two commuting sets of 
isospin generators: 

II dE4%i;o.i,E 

= J dEi/j^^Ei/^ii^E, (3.7) 

— * — * — # 

and similarly for I2. While the isospin generators Ii and I2, obey the usual SU{2) com- 
mutation relations, they commute with each other and also with the normal spin rotation 
generators, J (not to be confused with the coupling constants). 

In position space the one-dimensional problem can be defined in terms of left-movers 
only: 

Mt.^) = rf^e-^^(*+^/^Vi?, (3.8) 

^/V J~D 

for — 00 < X < 00. Here D is a bandwidth cut-off and v is an arbitrary velocity parameter 
which defines the scale of length in the effective one- dimensional problem. We adopt an 
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unconventional normalization for our one- dimensional position space fermion fields so that 
they obey: 

{V'i(a;),7/>z.(y)} = 27r5(x-y) (3.9) 

Alternatively, we may define right and left-movers on the half- line, x > 0, with 

,ljji{t,x)=^L{t,-x), (x>0) (3.10) 

and ipL{x,t) defined as in Eq. ( p.8|) . Note that Eq. ( |3.10|) would follow automatically from 
defining ipR^x^t) to be the right-moving field obeying the boundary condition ^/'^(t, 0) = 
iljiit,^)- (Henceforth we drop L, i?, subscripts. Fields are generally left-movers.) 

A key feature in the conformal field theory approach to this problem is bosonization. 
Actually, we use a somewhat generalized version of bosonization in which fermion fields are 
represented in terms of boson fields and Ising fields. A conventional abelian bosonization 
approach would involve introducing four bosons for the spin and charge degrees of freedom of 
each channel, (j)s,ii <Pc,i with i = 1,2. In fact we only introduce the left-moving components 
of these bosons, or equivalently introduce left and right movers on the half-line with an 
appropriate boundary condition. The bosonization formulas take the form: 

oc e^^('^=-"+<^='»), oc e^^('^-.'-<^=.») (3.11) 

As discussed elsewhere,^ these bosons should be regarded as periodic variables since only 
their exponentials (or derivatives) occur in physical quantities. Adopting a fixed normaliza- 
tion for the free boson Lagrangian density: 



C = U,ct>d^(t> (3.12) 



we define a "compactification radius" by identifying: 

<P = <P + 2'kR (3.13) 

The above bosonization formulas, Eq. ( p.ll| ), identify the compactification radius for the 
spin boson as i? = 1/\/2tx. At this radius, the free spin boson theory exhibits SU{2) 
symmetry and is, in fact equivalent to an SU{2) Wess-Zumino-Witten (WZW) non-linear a 
model with Kac- Moody central charge k = 1. The spin factors occurring in Eq. ( 3.11 ) can 
be written in terms of the left-moving factors of the WZW fields, Qaf- 



^Pa^^e'^'^'^-^ga^ (3.14) 

In fact, the charge bosons have the same compactification radius. This is no accident. 
It reflects the fact that the free fermion theories for each channel have an SU{2) x SU{2) 
symmetry. This corresponds to separate spin and isospin symmetry [see Eq. (p^)] for 
each channel. This is equivalent to 0(4). The 0(4) symmetry is manifest if the two spin 
components are each written in terms of hermitean and anti-hermitean parts, giving a total 
of four herimitean fermion fields, for each channel. The Kondo interactions will break the 
separate spin symmetries down to the diagonal subgroup; they have a similar effect on the 
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separate isospin symmetries except in the special case, J_ = 0. Thus we may introduce two 
more k = 1 WZW field, (/ij)^ corresponding to the two charge fields. Here A is an isospin 
index whereas i labels the two species. The bosonization formulas of Eq. then take 

the form: 

i^ai oc {hi)igai (3.15) 

The left-moving k = 1 WZW fields obey: 

^7"^ = e'^^gp (3.16) 

and similarly for the /ij's, as can be seen by comparing the abelian and non-abelian bosoniza- 
tion formulas. As we shall see, the charge bosons play a purely passive spectator role in 
the 2-impurity Kondo effect and it is not necessary to replace the charge bosons by the 
non-abelian /i-fields in what follows. 

The next crucial step is to rewrite the bosonized theory in terms of a total spin boson 
and some leftover degrees of freedom, describing relative spin fiuctuations. The obvious way 
of doing this, which is not the one we use, consists of changing variables to the sum and 
difference of the spin bosons, ips^i and 0s^2- The reason that we do not follow this procedure 
is that it does not explicitly maintain the total SU(2) symmetry of ordinary spin. Instead 
we use a procedure based on the Sugawara form of the Hamiltonian. 

Before explaining this in detail, we pause to review some basic properties of conformal 
field theories and Sugawara Hamiltonians. We work with left-movers only (so that all op- 
erators are functions of x + vt only) and scale a factor of v/2it out of the Hamiltonian for 
convenience, writing: 

H = ^l'd.H (3.17) 

Assuming that T-C{~1) = 'H{1)-, we define the Fourier transform: 



Ln = — dxe'^'^'^nix) (3.18) 



For a general conformally invariant Hamiltonian, these Fourier modes generate the Virasoro 
algebra: 

[L„, Lm] = (n - m)Ln+m + TTT^^^^ ~ l)'5n+m,0 (3.19) 

where c is the conformal anomaly parameter. One class of conformally invariant theories, 
relevant to the present discussion, has a Hamiltonian density quadratic in the currents, J°'{x) 
(a = 1, 2, 3, ...n) of some group of dimension n. In the present case, only the group SU{2), 
with n = 3, is relevant. The Fourier modes of the currents obey the Kac-Moody algebra: 

[J„^ J'J = te'^'^J^^^ + ^/cn<5„+„,o (3.20) 

Here k is the Kac-Moody central charge (or level). It must be a positive integer. The 
Hamiltonian then takes the Sugawara form: 
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^(^) = ^ ■ '^''(^) • ^''(^) ■ (3.21) 

or, in momentum space: 

oo 

Ln — 2 _|_ ^ ^ ^ • J—n ' Jn+m '■ (3.22) 

It follows from the Kac-Moody algebra that the Sugawara Hamiltonian obeys the Virasoro 
algebra with conformal anomaly, 

The Hamiltonian generates time-translation of the currents implying the additional commu- 
tation relations: 

[L„, rj = -mJ:^^ (3.24) 



This fixes the normalization of the Hamiltonian in Eq. (|3.21 



We now return to the problem at hand, considering at first the zero coupling, free theory. 
The free theory of two channels of left-moving fermions with spin can be written in Sugawara 
form. Altogether there are four commuting terms in the Hamiltonian, quadratic in the 
charge and spin currents for channels 1 and 2. The charge parts can actually be written 
more symmetrically as quadratic forms in the Isospin currents defined in Eq. ( p.6|) but this 
is not of great importance for what follows. The two sets of commuting spin currents for 
each channel obey the Kac-Moody algebra at level k = 1, and the associated Hamiltonians, 
Tisi, each have conformal anomaly c = 1. It is natural to rewrite the theory in terms of the 
total spin currents: 

J(x) = Ji(x) + J2(x) (3.25) 

These obey the Kac-Moody algebra with k = 2. The associated Hamiltonian, constructed 
from J using Eq. (|3.21|) with k = 2, Hg, has c = 3/2. A crucial point is that TCsi + 'Hs2 
and Tis both generate time-translation for the total spin currents. I.e., they both obey the 
commutation relations of Eq. ( |3.24|) with the total spin currents. Since Tig is itself quadratic 



in these total spin currents it follows that 

[Hsri^x) + ns2{x) -rLs{x),ns{y)] = o (3.26) 

Thus the spin Hamiltonian, Hsi + 'Hs2, can be written as a sum of two commuting pieces, TCs 
and a remainder. This is an example of the Goddard-Kent-Olive coset construction. § The 
remainder is associated with the coset [SU{2) x SU {2)]/SU{2) o in this case. ( "D" represents 
the diagonal subgroup.) The value of the conformal anomaly for the coset Hamiltonian is 
c = 1 + 1 — 3/2 = 1/2. The coset Hamiltonian obeys the Virasoro algebra with this value 
of c. There is a unique unitary conformal field theory with this value of c, namely the Ising 
model. 
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Thus rather than replacing the two spin bosons for each channel by a total spin boson 
and a difference boson we replace them by a /c = 2 Kac-Moody conformal field theory and 
an Ising model. The former theory can be considered to be the left-moving part of a /c = 2 
WZW model, ie., we may construct the various operators out of the left-moving factor of a 
unitary matrix, g^. The value of k for this matrix field is the coefficient of the Wess-Zumino 
topological term in the Lagrangian. We may now write down representations for the various 
operators in the free fermion theory as products of charge (or isospin) bosons, the total spin 
boson, Qa and the Ising field. The k = 2 WZW model has primary fieldfiJ of spin j = 
(identity operator, 1), j = 1/2 (fundamental field, go) and j = I (denoted (p.). The k = 1 
WZW model only has the identity operator and the j = 1/2 primary, h^. Their scaling 
dimension is given by the general formula: 

There are three primary fields in the Ising model: the identity operator, 1, the Ising order 
parameter, a of dimension x = 1/16 and the energy operator e with dimension x = 1/2. 
[These are dimensions of left-moving factors only. Thus the order parameter has total 
dimension 1/8 corresponding to a correlation exponent, rj = 2x = 1/4. The energy operator 
has total dimension x = 1 corresponding to a thermal exponent, u = 2 — x = 1.] The fermion 
field is written in this representation as: 

Ipai oc (hi)iga(T (3.28) 

Note that the three factors have dimensions which add up correctly to that of the fermion 
operator: x = 1/4 + 3/16 + 1/16 = 1/2. This determines uniquely the representation. Other 
operator representations can be determined using the operator product expansion (OPE). 
For the SU{2) WZW fields, hi and g, the OPE of two primary fields of spin j and j' gives 
each primary field with spin from \j — j'\ up to the minimum of j + j' and k — j — j'. For 
the Ising model the OPE gives: 

a X a ^ 1 + e, a x e ^ a, exe^l (3.29) 

This OPE is equivalent to that of the k = 2 WZW model with the identification of Ising 
and WZW primary fields: 

a g, e (p (3.30) 

Using the OPE, symmetry considerations and consistency of scaling dimensions, we can 
determine the representation of any operator in the free fermion theory. An important point 
is that not all products of isospin, spin and Ising operators occur in the representation of 
free fermions. When taking the OPE of a product of operators representing a free fermion 
operator, only a subset of all operators in the products of the OPE's for each factor occurs. 
This type of representation of free fermions is known as a conformal embedding. The 
complete set of products of primary fields that occurs in the representation of free fermions 
is given in Table I. Note that we only list primary fields. Other operators such as J are 
descendents. 
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A subtlety arises concerning parity. Parity takes ipo — > —ipo in the one-dimensional 
theory. In the alternative basis of Eq. (|3.3| ) it interchanges ipi and il)2- It thus follows that 
it interchanges the two isospin fields, hi and However, that is not the whole story as 
we see by considering the operator ipiail^i — ilj2^'ip2- This is odd under parity but has the 
representation 0e, independent of the isospin fields. To obtain a consistent definition of 
parity we define it to also take (p —(f). 

In the Bose-Ising representation the Kondo interaction is written as: 

Hk = J+/(0) ■ (5^1 + ^2) + J™0(O)e(O) ■ (5i - S2) + J_(/i1)^(O)(/i2)a(O)0(O) ■ (5i + ^2) 

(3.31) 

Here the three coupling constants, J±, Jm are proportional to the ones defined in Eq. ( ^.41) . 
Note that while J+ only couples the total spin field, Jm also couples the Ising field and 
J_ also couples the charge fields. Since J_ will turn out to be irrelevant for the particle- 
hole symmetric case being considered here, it follows that the relevant part of the Kondo 
interaction only involves the total spin and Ising fields. 

We may write the finite-size spectrum in terms of the "bosonized" representation. Here 
we take the one- dimensional effective theory and impose convenient boundary conditions on 
it on a line of finite length. Because of the rather unsymmetrical way that the dimensional 
reduction takes place [see Sec. II], these boundary conditions would not arise from any simple 
or natural ones on the original three-dimensional problem. Nevertheless we are interested 
in considering the finite system for two reasons. The first is that the Wilson numerical 
renormalization group method essentially studies a finite one- dimensional system. We will 
make detailed comparisons of the finite-size spectrum with the results of this method. The 
second reason is that there is an intimate connection between the finite-size spectrum and the 
operator content (for the infinite system) which we exploit. Thus we restrict the left-moving 
fermions to the interval: 

-l<x<l (3.32) 

and impose the boundary condition: 

V^l(0 = -^l(-0 (3.33) 

We note that in the formalism with left and right movers on the positive x-axis [see Eq. 
( p.lOp ] the boundary condition becomes: 

Ml)+Ml) = ^ (3.34) 

The free fermion spectrum is obtained by populating the single fermion momentum eigen- 
states with momentum k = {n + 1/2)tt/1. Note that with these boundary conditions, the 
groundstate is unique with all negative k states filled and positive k states empty, {k is 
measured from kp-) In general states have energies (measured from that of the ground- 
state) of the form E = {tt /l)x where x is an integer or half-integer. There is a one-to-one 
correspondence between states and operators with these boundary conditions; the value of 
X corresponding to the dimension of the operator. For instance, x = 1/2 corresponds to 
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a single particle or hole state. It is created from the groundstate by a single application 
of the fermion field. Just as the operators can be written as products of charge, spin and 
Ising operators, the states can be written as direct products of states from each sector. In 
each sector, the states are grouped into several conformal towers with energies of the form 
E = {'k/1){x + n) where x is fixed and fractional and n is an integer. The conformal towers 
are in one-to-one correspondence with the set of primary fields, with x being the scaling 
dimension. The set of products of conformal towers from each sector which can occur in 
the free fermion spectrum corresponds to the set of products of primary fields that occur in 
representing free fermion operators. Thus it can be read off from Table I. 

IV. FINITE-SIZE SPECTRUM 

In this section we postulate a critical theory of the unstable multicritical point which 
occurs in the particle-hole symmetric case at a critical value of the inter-impurity coupling, 
K. We derive the corresponding finite-size spectrum and compare it to the results of the 
NRG. 

A fundamental assumption behind our approach is that the critical point can be described 
as a conformally invariant boundary condition.0 We expect this to be generally true for a 
wide class of critical phenomena involving quantum impurities. Indeed this assumption is 
very analogous to the one which is widely made in studying bulk critical phenomena in 
two-dimensional [or (l-l-l)-dimensional] systems. There it is assumed that a wide class of 
critical phenomena exhibits conformal invariance. Although such an assumption can rarely 
be proven, it is consistent in the following sense. 

We first of all assume scale invariance at the critical point. This could be taken as a 
definition of criticality. We then argue that the critical theory should have S0{2) or Lorentz 
[SL{1, 1)] invariance. Although this is not a symmetry of most underlying microscopic 
theories it can be seen that the operators which break this symmetry down to whatever 
subgroup exists in the microscopic theory are irrelevant. In the case of a two-dimensional 
theory defined on the square lattice, the subgroup of 5*0 (2) is the symmetry group of 
the square lattice. This permits non-rotationally invariant derivative terms in a Landau- 
Ginsburg theory such as {dcp/dx)^ + {d(j)/dy)^. All such terms are irrelevant. For a one- 
dimensional quantum fermion system the dispersion relation is not usually exactly linear. 
However, the additional terms in the Hamiltonian reflecting the non-linearity involve at 
least two derivatives and hence are usually irrelevant. (Of course the situation changes if 
the linear term vanishes.) Once we have convinced ourselves that our critical theory is scale 
invariant and Lorentz invariant it follows immediately that it must be conformally invariant. 

A quantum impurity problem cannot be invariant under the full conformal group. Indeed 
it is not even Lorentz invariant since a special point is singled out; or a special line, x = in 
the space-time description. The maximal symmetry that we could hope for is the subgroup 
of the conformal group which leaves this line fixed. In the imaginary time formulation, 
introducing the complex co-ordinate, z = vr + ix, the full conformal group is all analytic 
transformations z ^ w{z). The subgroup leaving the line a; = 0, (the real axis) invariant is 
the set of transformations for which w{T)eTZ. Taylor expanding the analytic function: 
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oo 
n=l 



(4.1) 



we see that while, for a general conformal transformation the a„'s can be arbitrary complex 
numbers, in the presence of a boundary they must all be real. Thus, the (infinite) number of 
symmetry generators is reduced by a factor of 2. This infinite set of symmetries corresponds 
to time translations and rescalings of space and time which may be performed independently 
at each space-time point. 

We generally expect that a system which exhibits conformal invariance far from a bound- 
ary will exhibit boundary conformal invariance as the boundary is approached. Such a sys- 
tem could be a two-dimensional statistical system at its bulk critical temperature with a 
boundary. It could also be a one- dimensional quantum system in a gapless phase with linear 
dispersion relation in the presence of a quantum impurity. Note that we have managed 
to formulate the two-impurity Kondo model as such a system. In the left-right formalism 
the model is defined on the positive x-axis with both impurities at x = 0. Of course, the 
universal critical behavior only emerges in the scaling limit. If we calculate correlation func- 
tions we only expect universal behavior when all points are well separated. Furthermore all 
points must be far from the boundary compared to microscopic scales. Since the ratios of 
distances from the boundary to distances between the points remain as free parameters it 
is still possible to observe critical behavior which is affected by the boundary. Generally we 
can also make universal predictions when the points are close to the boundary but far from 
each other. However in this case only exponents, not amplitudes are universal. Again if we 
assume invariance under time translations and global rescaling of space or time the local 
invariance follows. 

The specification of conformally invariant boundary conditions can be conveniently for- 
mulated in terms of the finite-size spectrum. We consider a system defined on a cylinder of 
length / and circumference /3 (ie. a quantum system at temperature T = 1//3) with some 
conformally invariant boundary conditions A and B at the two ends. Conformal invariance 
implies that the spectrum can be specified in terms of the conformal towers of the periodic 
left-moving system on an interval of length 21. The partition function can only be a sum of 
partition functions for each conformal tower with a multiplicity factor: 

Z{l/f3)=Y.n\sm/P) (4.2) 

i 

Various pairs of conformally invariant boundary conditions correspond to various sets of 
integers n^^. It turns out that not all possible choices of integers correspond to a pair 
of conformally invariant boundary conditions. Cardy derived a set of powerful consistency 
equations that these integers must obey.lii In the particular case where the theory is a 
product of several decoupled theories, as for all versions of the Kondo problem, the set of 
conformal towers in Eq. (|4.2| ) must be summed over all products of conformal towers from 
each sector. For the two-impurity Kondo problem, there are four sectors: the two isospins, 
total spin and Ising. The number of conformal towers in each sector is two for each isospin 
(labelled by the isospin of the primary field, i = or i = 1/2) three for spin {j = 0, 1/2, 1) 
and three for Ising (1, a and e). Thus altogether there are 36 products of conformal towers 
and 36 integers to specify. For the free fermion spectrum discussed in the previous section, 
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six of these integers have the value 1 and the rest are 0. The 6 products of conformal towers 
that occur are given in Table I. 

It turns out that a useful way of generating new conformally invariant spectra from old 
ones is "fusion" This corresponds to a particular mapping of each conformal tower into a 
set of other ones. The fusion rules are not obviously related to boundary critical phenomena 
per se but come from the bulk operator product expansion. If we consider the OPE of two 
bulk primary operators, Oi and Oj then it will contain the complete set of primary fields 
Ok with multiplicities N^j. The new spectrum is obtained by the replacement: 

^'-E^>' (4-3) 

k 

for any fixed j. This procedure gives a new solution of Cardy's consistency equations for any 
choice of the conformal tower, j. The fusion rule coefficients are a property of each sector 
of the theory independently. They were given for the Ising sector in Eq. ( 3.29| ). 



As a trivial example of a new conformally invariant spectrum generated by fusion we 
consider fusion with the j = 1 field in the spin sector, 0. We see from the fusion rules that 
this simply interchanges the j = and j = 1 conformal towers. The number of products of 
conformal towers remains six and is given in Table II. Notice that in this case the groundstate 
is 16-fold degenerate. This spectrum corresponds to free fermions with a 7r/2 phase shift.0 
I.e. the Fermi level now sits at one of the single-particle energy levels. Since each of the 
four species of Fermions (two channel times two spins) can have this level filled or empty, 
this accounts for the 16-fold degeneracy. In fact precisely this spectrum describes the Kondo 
fixed point which occurs for large ferromagnetic inter-impurity interaction. We also note, 
that the same spectrum is obtained by fusion with either the e field in the Ising sector or 
the ii = i2 = 1/2 field in the isospin sector. 

The other possible fusion in the Ising sector is with the a field. This gives the non- 
trivial critical point in the two-impurity Kondo problem. The resulting spectrum is given 
in Table III. In this case there are eight products of conformal towers. Note that the same 
spectrum can be obtained by beginning with the spectrum after 0-fusion, of Table II and 
then performing a-fusion. Thus the non-trivial fixed point is symmetric with respect to 
the zero phase shift and 71/ 2 phase shift fixed points. We now wish to compare this in 
detail with the NRG results. This is done by multiplying together the various conformal 
towers. A general state will have energy E = {'kv/1){x + n). Here x is the energy of the 
primary (or highest weight) state and the integer n is the descendent level. Both x and n 
are obtained by summing over the four sectors. The values of x that occur for the eight 
products of conformal towers are given in Table III. We will measure energies from that of the 
groundstate; thus x — 1/16 occurs. We will content ourselves with working out the energies 
of all states with x + n — < 2. This implies n = or 1 only. Thus we only need know 
the first descendents in each conformal tower. Altogether there are 8 different conformal 
towers in the three inequivalent component theories {k = 1 WZW, k = 2 WZW and Ising). 
All first descendents of all 8 conformal towers are given in Table IV. We now obtain the 
complete spectrum at the non-trivial fixed point with x + n — l/16<2by taking either the 
primary state in each sector or taking a first descendent in one sector and primaries in the 
other three. This gives the set of states shown in Table V. Recall that parity interchanges 
the two isospins and also multiplies by (—1) for the conformal tower. The spectrum is 
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symmetric between the two isospins. By taking symmetric or antisymmetric products of the 
two isospins we obtain states of definite parity. 

The general particle-hole symmetric Hamiltonian ((2.6) or (2.9), with arbitrary coeffi- 
cients) has three symmetries: total isospin / = Ii + I2, total spin, and (exchange) parity. 
When the coupling constant J_ = 0, a further symmetry develops, and one can describe the 
eigenstates with an alternative set of three quantum numbers: the two isospins, ii and i2, 
and total spin. This was actually the case in the NRG calculation to which we compare the 
conformal field theory (CFT) spectrum. Eigenstates within the same ii, 12 isospin-parity 
multiplets are seen in Table V to be degenerate to better than one part in 10^, an indication 
of the accuracy to which the symmetry holds in the NRG calculations.^ 

In Table V we express the lowest eigenstates via both sets of quantum numbers. Columns 
one and two constitute a complete set: the two isospins [iiM) and the total spin (j). To 
facilitate comparison with the NRG, we have also resolved the spectrum into states of 
definite i (total isospin), j (total spin), and P, parity. (Note that resolving states of definite 
i picks out states of definite symmetry with respect to interchanging the two isospins and 
hence determines states of definite parity. In particular the z = part of an (^1,^2) = 
(1/2, 1/2) multiplet is antisymmetric and the i = 1 part is symmetric. An additional intrinsic 
odd parity must be included for all states whose spin factor is in the j = 1 conformal 
tower.) Column five of Table V indicates i and P (+ = even, — = odd) in this alternative 
formulation. Altogether we obtain 49 multiplets with definite {i,j,P). In the last column 
in Table V we give the energies of the lowest states of the same quantum numbers obtained 
using the NRG. 

Very briefly, numerical renormalization group calculations" iteratively diagonalize a 
Hamiltonian, bringing in successive low-energy degrees of freedom until asymptotically the 
ground state is reached. After each iteration, the set of energy levels is ordered and trun- 
cated to keep only a fixed number of the lowest-energy states. Truncation errors are the 
main source of inaccuracy in the energy levels, and although the truncation approximation 
is expected to be a good one for the very lowest energies, the uncertainty (typically a few 
percent at most) increases with increasing energy. For the results in column six of Table 
V, approximately 1100 states were kept at each iteration. The coupling constant uJ = 0.18 
(c.f. Eqs. (2.6, 2.7)). For this value of J, the unstable fixed point is fully realized by 14 
iterations. 

To determine the critical point, the ratio of RKKY coupling to Kondo temperature was 
varied until the region of unchanging energy levels was maximized. In practice, of course, 
numerically one can never exactly sit at the critical point; since the critical point occurs for a 
finite (probably irrational) value of the coupling parameters, it is in practice only possible to 
choose initial parameters such that the flows are asymptotically close to the critical point. 
And, in fact, no matter how carefully one adjusts the parameters, at any iteration there 
are always some energy levels that are not completely fiat. Since for an inexact choice of 
parameters the critical region does not extend all the way to T = 0, one must choose a 
temperature (iteration number) at which the approach to the critical point is closest. These 
two facts, that one is never right at the unstable point in parameter space, and that even 
with a very good initial parameter choice, it is not clear which iteration to choose to best 
represent the unstable fixed point, thus introduce additional uncertainty in the numerical 
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energy values in Table V. 

The agreement with NRG nonetheless seems very satisfactory. The 49 NRG states shown 
in Table V consist of numbers 1—47, 68 — 69 in a consecutive numbering by increasing energy. 
The overall energy scale, set by f , is not determined by the NRG so we have adjusted it 
to fit the CFT spectrum. In a few cases the CFT predicts several multiplets with identical 
quantum numbers and the same energy, so the matching with the NRG levels involves some 
arbitrary choices. Otherwise, there are no free parameters. The agreement is to within 
about 2% for the first 17 levels but gets progressively worse at higher energies. At the 
highest energies shown the disagreement is about 5% for most of the levels. @ The two sets 
of multiplets involving the | descendent in the spin sector (column 2 in table V) have 
unusually large disagreements with CFT compared to other states at the same energy, 11% 
and 19% respectively. It may be possible to understand this in terms of the contribution of 
the leading irrelevant operator to finite-size effects.0 We regard Table V as rather convincing 
evidence that we have identified the correct boundary CFT critical point for this problem 
and that the NRG method is highly accurate. 

It seems clear that the spectrum of Table V could not be obtained from simply imposing 
linear boundary conditions, respecting spin-rotation symmetry, on the free fermion fields, 
and filling up the corresponding levels using Fermi statistics. Note, for example, the strange 
ratios of energy gaps: 3:4:7:8 and the peculiar 7-fold degeneracy of the states with 
X — 1/16 = 1/2. Thus we refer to this as a "non- Fermi liquid spectrum". Nonetheless, 
it is possible to describe this spectrum in terms of free Major ana Fermion field and their 
spinor representations, using the alternative Ising xS0{7)i bosonization scheme of Ref. ( 



141) . These fields are related in a non-local way to the original conduction electron degrees 



of freedom. A related physical picture can also be found in Ref.'s ( |T5|,|T6D. We explain this 
hidden symmetry in detail in Sec. V. 



V. OPERATOR CONTENT 

The operator content at the non-trivial fixed point is determined by "double fusion" 
with a. This is because we must effectively consider the spectrum on a strip with the 
same boundary condition at each side. The result in shown in Table VI. To understand the 
stability of the non-trivial fixed point we must carefully consider all relevant and marginal 
boundary operators {x < 1) which are allowed by symmetry to appear in the Hamiltonian. 
(As usual in critical phenomena we expect that "anything that can happen will happen" .) 
Apart from the identity operator, there are seven other multiplets of relevant operators with 
X = 1/2. If we assume that J, Ji and I2 are all conserved then only the seventh operator 
in Table VI, which we may write simply as e, may appear in the Hamiltonian. The J_ 
coupling constant of Eq. ( p. 31 ) breaks isospin down to the diagonal subgroup. This might 



seem to permit an element of the eighth operator multiplet in Table VI: {h\)'^{h2)A- Using 
Eq. ( |3.16|) , this can be rewritten as (hi) A{h2) b^^^ ■ Since we may regard the fields hi and h2 



as commuting, it follows that this operator is odd under parity. This conclusion also follows 
from the form of the Hamiltonian in Eq. ( p.31|) ; since (j) is odd under parity, so must be 



{h\)^{h2)A- Similarly the various marginal primary operators in Table VI are not allowed 
by symmetry. There are also marginal descendent operators, the currents J and Jj. These 
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are forbidden by spin and (diagonal) isospin symmetries. 

Hence we reach the important conclusion that there is only one relevant operator permit- 
ted, e. It follows that the non-trivial critical point can be reached, in general, by tuning one 
parameter to its critical value. We could vary any of the coupling constants in the problem 
to reach the critical point. The physical picture is clearest if we vary the RKKY coupling 

as explained in the introduction. We expect that varying K away from its critical value, 
Kc, drives the system to one of the stable fixed points with either a 7r/2 phase shift in both 
channels or zero phase shift in both channels. The coefficient of e in the effective Hamilto- 
niam should be proportional to A' — K^.. In particular, it is important to note that for the 
case of particle- hole symmetry, tuning such that J_ = 0, that is, setting even and odd-parity 
couplings equal, is not necessary for obtaining the critical point as long as K = Kc. 

There is a formal analogy between this renormalization group flow and one that occurs 
in the classical Ising model.EJ Consider a two dimensional Ising model at its critical temper- 
ature, defined on a half-plane with free boundary conditions. Then apply a weak magnetic 
field at the boundary only, pointing up or down. The free boundary condition represents an 
unstable fixed point; applying the magnetic field drives the system to the stable boundary 
condition with spins pointing up or down at the boundary. The local magnetic field corre- 
sponds to the e operator of dimension 1/2. The three different boundary conditions, spin 
up, spin down and free are related to each other by fusion, e fusion takes spin up into spin 
down and vice versa; a fusion takes spin up or down into free. 

We may approach the critical point in various ways; one is by setting T = and letting 
K approach K^, another is by setting K = Kc and letting T go to zero. If we tune K to Kc, 
then the corrections to scaling behaviour are governed by the leading irrelevant operator. 
None of the other primary operators in Table VI are allowed by symmetry. The leading 
irrelevant operator is the x = 3/2 first descendent e'. It can be checked that the various 
other dimension 3/2 operators are not allowed by symmetry. In particular, J_i ■ 0, the 
leading irrelevant operator in the multi-channel, single-impurity Kondo effect is not allowed 
by parity. 



We may define the coefficient of the e' term in the effective Hamiltonian as 1/ \/Tk, where 

Tk is the Kondo temperature. Thus, at low temperatures and for K close to Kc, we must 
consider the following Hamiltonian: 



where a is a constant and Hpp is the Hamiltonian at the fixed point, with the non-trivial 
boundary condition. We first consider the properties of Hpp itself, assuming that K = Kc 
and ignoring the effects of the irrelevant operator. We then consider a non-zero e' with e still 
zero, corresponding to approaching the critical point by setting K = Kc and T ^ 0. Finally 
we consider e and e' both non-zero corresponding to a general approach to the criticahpoint. 

A quantity which characterizes the critical point is the residual impurity entropy.EJ This 
is defined by taking the length of the system to oo first and then taking the temperature 
to zero. In general the impurity entropy, i.e. the part which does not scale with length, 
but rather goes to a constant as the length goes to oo, goes to a finite non-zero constant 
as T ^ 0. This quantity is always the log of an integer for a Fermi liquid fixed point, the 




H = Hpp + a 




e + 




1 



(5.1) 
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integer simply being the degeneracy of the decoupled impurity. In the two-impurity Kondo 
problem, we expect zero degeneracy at the stable fixed points since the impurity spins are 
always in singlet states. At a non- Fermi liquid fixed point, the residual impurity entropy, 
S{0) = In^f where g need not be an integer, in general. When the fixed point is obtained by 
fusion, a general expression for g has been derived in terms of the modular S-matrix. This 
gives: 

g = ^^=V2. (5.2) 

Note that g decreases under renormalization from the unstable critical point to the stable 
ones, in accord with the "(/-theorem" . 

Another striking signature of the non-trivial critical point is provided by the zero-energy 
one-particle S-matrix, 5(1) .1 This is the amplitude for a single electron to scatter off the 
impurity and produce a single electron final state. If the system exhibited Fermi liquid 
behavior then S'(i) would have unit amplitude since all multi-particle final states would have 
vanishing probability at zero energy. Thus S'(i) would reduce to a phase shift, e^*''. However, 
at a non- Fermi liquid fixed point, the probability of inelastic scattering (one particle into 
many) does not vanish at zero energy so S'(i) can have non-unit amplitude without violating 
unitarity. 1 — provides a direct measurement of how "non- Fermi liquid like" the critical 
point is. It is a universal quantity characterizing the various fixed points. It is defined 
precisely by the asympotic behaviour of the single electron Green's function. In the absence 
of any impurity couplings, this would have the non-interacting form: 

<^f{z)^,,{z')>^^ (5.3) 

This dependence on spin indices will always be of the above form, due to spin-rotation 
symmetry, so we will suppress it in what follows. Ignoring the irrelevant, J_ term in Eq. 
( p.4| ), the same is true of the dependence on the channel indices. If x and x' have the same 
sign, the form of Eq. (|5.3|) holds even in the presence of impurity interactions. This is a 
general feature of conformally invariant boundary conditions. On the other hand, when the 
two fields at x and x' straddle the impurity, then different asympotic behaviour can emerge, 
which depends on the universality class of the boundary interactions. The most general 
form consistent with conformal invariance and spin and channel symmetries is simply: 

<Cn-)V^..(-')>-^^ (5.4) 

where S'(i) is a universal amplitude, which we can interpret as the S-matrix element for one- 
particle into one-particle scattering. In the non-interacting case, 5(1) = 1. This should also 
be the case at the stable Fermi liquid fixed point where the two impurities form a singlet. On 
the other hand, at the Fermi liquid fixed point where both impurities are Kondo screened, 
^(i) = —1, corresponding to a 7i/2 phase shift. As mentioned in Sec. I, 5(1) must be real in 
all cases by particle-hole symmetry. 5(i) is given by the matrix element of the a operator 
(since this is the Ising factor of the fermion operator, ip) between the boundary state and 
the vacuum state. Since the boundary state is obtained by fusion with the a primary field, 
this gives the result: 
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= (5-5) 



where Sj is the modular S-matrix (which has no obvious connections with the scattering 
matrix) of the Ising modeL It turns out that = 0, while the other matrix elements are 
non-zero. Hence, we conclude that: 

5(1) = (5.6) 

at the non-trivial fixed point. This implies that the scattering is entirely inelastic at zero 
energy. Thus, is a sense, this fixed point is as non-Fermi liquid like as possible. The fact 
that 5(1) = can be understood as a consequence of the symmetry between the trivial fixed 
points with 5(1) = 1 and 5(i) = — 1. As mentioned in Sec. Ill, the non-trivial fixed point 
can be reached by cx-fusion starting from ezt/ier trivial fixed point. This implies that 5(i) at 
the non-trivial fixed point is proportional to its value at both trivial fixed points with the 
same constant of proportionality. Since 5(i) has opposite sign at the two trivial fixed points, 
it follows that it must vanish at the non-trivial one. 

If we relax the assumption that J_ = 0, then the Green's function need no longer be 
diagonal in channel indices. We may also relax the assumption of particle-hole symmetry. 
In general, at a Fermi liquid fixed point we may write: 

<^l.f{z)^,,iz')>^^-^, (5.7) 

where 5(1),^ is a general unitary matrix. It is the non-unitarity of 5(i) which signals non- 
Fermi liquid behaviour. 

Ignoring all (relevant or irrelevant) perturbations about the critical point, the impurity 
specific heat and the uniform impurity susceptibility, which is the response of the system to 
a conserved quantity, vanish. This occurs for the same reason as in the single impurity case 



and we refer the reader to Ref. ( |TT]). On the other hand the staggered susceptibility does 
not vanish. This is defined in terms of the correlation function of 5i — This operator is 
odd under parity so we expect that it reduces to the primary field at the critical point. 
This operator has dimension 1/2 giving a correlation function: 

< [Si - S2] (r) • [Si - S2] (0) >oc ^ (5.8) 

|r| 

This is the same as the impurity spin correlation function in the overscreened two-channel 
single impurity case mentioned above. Fourier transforming at finite temperature and an- 
alytically continuing to real frequency, we obtain [see third reference in [ the staggered 
susceptibility. The real part behaves as: 

Rexs(^, T) oc ln(T^/T) {u « T) (5.9) 



Tk 
1 

Tk 



^^\n{TK/uj) {T«uj). (5.10) 



The imaginary part behaves as: 
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Imxs oc tanh(/3u;/2). (5.11) 

The electron spin operator with the same symmetry, iplaipi — ip2^ip2 should also couple to 
and exhibit the same decay, close to the boundary. The universal crossover function from 
bulk to boundary behaviour could be calculated following the techniques of Ref. ( |^). Various 
other operators exhibit the same |r|~^ decay due to couphng to other x = 1/2 boundary 
operators in Table VI, e and {hi)A{h2)B- As already remarked, the various interaction terms 
in the Hamiltonian such as Si ■ §2 couple to e and so should have this decay. Similarly, 
the spin-singlet, channel symmetric pair operator, {i^ai4'i32 + '?/'a2'?/'/32)e"^ should couple to 
(/ii)i(Mi. 

Now we consider the effect of the leading irrelevant operator, (l/-\/7^)e', still assuming 
that K = Kc- We can calculate physical properties by doing perturbation theory in l/y/Tx- 
This is similar to the calculations for the two-channel single impurity, s = 1/2 Kondo effect, 
where the leading irrelevant operator also has dimension 3/2. However, in this case the 
leading irrelevant operator is a Virasoro descendent, rather than a primary. Indeed, we may 
consider it to be the derivative of the primary field, e: e' ~ de/dr. Consequently the finite 
temperature two-point function has a different form than in the two-channel case and we 
find that there is no logarithmic singularity in the specific heat. I.e., 

^imp °^ T/Tk- (5.12) 

Note that the uniform susceptibility [ie. the response to a field coupling to the conserved 
total spin of impurities plus conduction electrons] obtains no contribution to second order in 
perturbation theory in the e' term and therefore also has no logarithmic singularity. This is 
again unlike the two-channel case where another dimension 3/2 irrelevant operator, J_i ■ 
is present in the effective Hamiltonian. It is absent in the two-impurity case by parity. We 
expect a finite impurity susceptibility as T — 0, of 0(1/Tx), coming from the dimension 
two irrelevant operator J^. 

Let us now consider the behavior for K slightly displaced from Kf.. As mentioned, we 
expect the system to renormalize to one or the other of the stable Fermi liquid fixed points. 
At these fixed points the impurity spins either form an inter-impurity singlet or else are 
Kondo-screened. There are various leading irrelevant operators of dimension 2. One such 
operator comes from the Ising sector; it is the Ising energy- momentum tensor, T/. When K 
is very close to Kc the system does not renormalize away from the unstable critical point 
until a very low energy scale is reached of order [K — K^Y /Tk- This plays the role of the 
lowest effective Kondo temperature at the stable fixed point. Thus we expect a term of the 
form 

in the effective Hamiltonian at the stable fixed point. The other leading irrelevant operators, 
such as in the spin sector, presumably have coefficients which remain finite as K ^ K^. 
We expect this because the delayed flow away from the unstable critical point only occurs 
in the Ising sector. The uniform impurity susceptibility goes to a constant as T ^ and 
can be determined by first order perturbation theory in the leading irrelevant operators. It 
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should remain finite as K ^ Kc, because it doesn't involve any Ising sector fields and is 
determined by the coefficient of J^. On the other hand, the specific heat does get a first 
order contribution from the irrelevant operator of Eq. ( p.l3| ). Therefore, we expect a specific 
heat slope. 

In general, for finite T and finite K — Kc, we expect the specific heat to be given by some 
universal scaling function: 

_ T /T K-K, \ 

where the scaling function f{x,y) has the asymptotic behaviour: 

/(x, 0) — » constant (x 0) 

constant , , , , 

f{0,y)^ ^ y-O. 5.16 

These scaling results are consistent with those of the NRG, which, at T = finJ that 
Cimp/T diverges as 1/{K — Kc)^, and that the uniform susceptibility is non-singular as 
K ^ Kc, in the presence of particle-hole symmetry. 

We now wish to consider the effect of particle-hole symmetry breaking. The first type of 
P-H symmetry, defined in Sec. II corresponds to: 

tpai eaf3^f^ (5.17) 

In the bosonized representation, this has the effect (^i)^^ or equivalently: 

{hi)A eAB{hi)B (5.18) 

The relevant particle-hole symmetry breaking interactions correspond to potential scattering 
terms of the form: 

SH = Ve4ije + Voi^l^o (5.19) 

(Here, the operators are evaluated at the origin in position space.) In terms of the fields, 
ipi, this becomes: 

SH = ^^^^[4^1 + 4h + ~ [4^2 + 4H (5.20) 

Upon bosonizing, this gives: 

6H oc (K + Vo){I' + I2) + constant ■ (K - K)tr(/ii)V3(/i2)e (5.21) 

The first term produces the marginal current operator. In terms of an abelian total charge 
boson, (j), this is proportional to d(j)/dx. We expect this to be exactly marginal and to 
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produce a line of fixed points with continuously varying phase shifts connected to the non- 
trivial critical point. This is similar to adding potential scattering to the two-channel s = 1/2 
critical point. We expect the second operator to produce the relevant operator tr(/ii)V^(/i2) 
at the non-trivial critical point. There are various ways of seeing that this should occur. One 
is by observing that, after double fusion, e can turn into the identity operator. Another is 
simply to observe that this relevant operator transforms the same way as the second term in 
Eq. ( [5. 211 ) under all the symmetries of the problem. IfVe = Vo and J_ = 0, then the relevant 
operator is not allowed by symmetry, since the original Hamiltonian has two independent 
symmetries If and /| whereas the relevant term only respects the diagonal symmetry If + I^- 
If J„ 7^ then the relevant term is allowed by symmetry even if Vf. = Vq. In this case we 
expect it to appear. Thus we conclude that potential scattering is only relevant if Ve Vo 
or J_ 7^ 0. We expect that the relevant potential scattering term will drive the system to 
a Fermi liquid fixed point, corresponding to a unitary S'(i) of the general form of Eq. ( ^.Tj) . 
In the limit where the relevant potential scattering term, tr(/ii)"'"r'^(/i2) has a very small 
coefficient, the flow to the stable, Fermi liquid, fixed point occurs at a very low energy scale, 
so we again expect that one of the leading irrelevant operators at the stable fixed point 
will have a large coefficient of 0[TK/{Ve — Vo)"^] (assuming that J_ =0). This time, the 
large coupling constant should occur in the isospin sector, rather than the Ising sector, as 
above. Thus the relevant potential scattering term should produce a diverging specific heat 
coefficient oc 1/(K — VoY- We do not expect uniform or staggered spin susceptibilities to 
diverge as K — Vq — > in this case, since the spin operators do not involve the isospin fields. 
However, divergences should occur in charge and pair correlation functions. 

Earlier NRG calculationgH^ are all in agreement with these predictions of CFT regarding 
particle-hole symmetry breaking. When J_ = and Vg = V^, the two Fermi liquid fixed 
points and the non-Fermi liquid critical non-trivial critical point are each extended into a line 
a fixed points with continuously varying phase shift. Potential scattering is thus irrelevant 
in this case. When J_ 7^ or Ve 7^ Vo, all T = flows are to a single line of Fermi liquid 
fixed points, connecting the antiferromagnetic RKKY singlet at 5 = to the two-impurity 
Kondo point at 5 = 7r/2. Potential scattering is relevant, and the non-trivial critical point 
is washed out in this case. 

Note that the first type of P-H symmetry is necessary to stabilize the non-trivial fixed 
point. The second type allows potential scattering with Ve = — Vo- This permits the relevant 
interaction discussed above [the second term in Eq. ( ^.21|) ] so the non-trivial critical point 
should not occur in this case. 

Finally, we show that an anisotropy of the Kondo exchange coupling constants J-t, 
gives rise only to an irrelevant perturbation of the non- Fermi-liquid fixed point. This shows 
that the exchange-anisotropic version of the 2-impurity Kondo Hamiltonian gives rise to 
the same universal low temperature physics as the isotropic one, described above. More 
precisely, we consider an anisotropic Kondo exchange coupling, such that the Hamiltonian 
is only invariant under rotations about the z-axis in S'f/(2)-spin space (so far, it has been 
invariant under rotations about any such rotation axis). To show the irrelevance of the 
exchange anisotropy, we proceed as in Ref.( |l^). Since the exchange anisotropy breaks the 
total spin SU{2) symmetry, some operators may occur in the fixed point hamiltonian which 
were not allowed in the isotropic case on symmetry grounds. We will now show that the 
only additional such operators are irrelevant. To see this we simply need to consult Table 
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VI, of boundary operators which exist at the non-Fermi-hquid fixed point: Operators of 
spin j = 1/2 are not allowed since rotations about 27r about any axis in SU{2)-spin space 
leaves the anisotropic hamiltonian invariant, but multiplies these operators by (—1). Thus 
we are left with the spin j = 1 operators in Table VI: Amongst those, the only relevant 
or marginal new operators (scaling dimension a; < 1) that could occur in the anisotropic 
case are (i): (f)^, with quantum numbers {ii,i2, j, Is) = (0,0,1,1), (ii): 0^e, with quantum 
numbers (0, 0, 1, e), as well as (iii) (/ii)'^(/i2)A0^ with quantum numbers (1/2, 1/2, 1, 1) and 
(iv): the spin-descendant J^, with quantum numbers (0,0,1', 1). However, operators (i) 
and (ii) are odd under parity, and are thus not not allowed to occur since the microscopic 
(exchange anisotropic) Hamiltonian is parity invariant. Operator (iii) changes sign under 
rotations by vr about the x— or y— axis in spin space whereas the microscopic Hamiltonian 
does not. Hence (iii) cannot occur. Operator (iv) is odd under time reversal [this can be 
seen e.g. from Eg. (|3.31| )] and thus cannot occur either ( the latter does, in any case, not 
renormalize, and could not destroy the non- Fermi-liquid fixed point for that reason). Hence 
we conclude that no relevant operators can occur at the fixed point even in the presence of 
exchange anisotropy. 



VI. HIDDEN SO (7) SYMMETRY 

An interesting feature of the spectrum at the non-trivial fixed point, discussed in Sec. Ill, 
is various apparently "accidental" degeneracies. For instance the second and third conformal 
tower in Table III have the same value of x — 1/16 = 3/8 although they have different values 
of {ii,i2,j)- This suggests the presence of some type of symmetry which is higher than the 
spin and isospin symmetries discussed so far. A related observation was made in Ref. ( p4D : 
there, an additional conserved quantity was noticed at the non-trivial fixed point, namely, 
/ dx [Je — Jo^ where, 

X- Jo = ipl^^e - ^l^^o = V'i|V'2 + V'^l^i- (6.1) 

In fact, a higher symmetry has been found in Ref. ( 0) and corresponds to the group S0{7) 
which contains SU{2y^ x SU{2Y^ x SU{2y as a subgroup. [ This is in fact a conformal 
embedding, S0{7)i SU{2){' x SU{2)\' x SU{2)i.] In this section we spell out this 
symmetry in more detail. 

We emphasize at the outset that this symmetry is rather different and more mysterious 
than the isospin and spin symmetries discussed so far in that it does not become manifest for 



any choice of the Kondo couplings in Eq. ( |3.4| ), except the trivial case where all couplings 



are zero. In this case the full symmetry is actually 5*0(8). This follows because we have 
4 species of complex or Dirac fermions, equivalent to 8 Hermitean or Majorana fermions. 
This 50(8) symmetry is also present at the Kondo screened fixed point (where the RKKY 
interaction is strongly ferromagnetic) since the phase shift boundary condition does not 
break the symmetry between the 8 Majorana fermions. Remarkably, an 5*0(7) subgroup 
of this 5*0(8) also remains at the non-trivial fixed point. The key to understanding this 
symmetry is the realization that the non-trivial critical point is related to the free fermion 
one by fusion in the Ising sector. Thus what we must do is to represent the free fermions as a 
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product of an Ising model and another sector of maximal possible symmetry. This symmetry 
will be preserved under fusion in the Ising sector. The maximally symmetric representation 
corresponds to a conformal embedding or alternative non-abelian bosonization scheme where 
the 8 free fermions are represented by an S0{7)i Kac-Moody conformal field theory together 
with an Ising conformal field theory. This preserves the value of the conformal anomaly or 
Virasoro central charge parameter c, since the KM algebra 5*0(7)1 has c = 7/2 and the Ising 
model, c = 1/2: 

c = 4 = 7/2 + 1/2. (6.2) 

There are actually two inequivalent ways of constructing a conformal embedding depending 
on the SO (7) representations into which the 8 Majorana fermions, transforming under the 
vector representation of SO{S), decompose. The trivial embedding corresponds to: 

S0{8)i S0{7)i X Is 

8^(7,l) + (l,e). (6.3) 

[Here and in the following we label SO{n) representations by their dimensions and conformal 
towers by the representation of their highest weight state.] This simply takes advantage of 
the fact that the Ising model is itself equivalent to a free Majorana fermion. We may 

fh 

define the 5*0(7) transformation to rotate 7 of the fermions and leave the 8 invariant. 
However, this is not the appropriate embedding; ie. this Ising model does not appear in the 
representation of the fermions in the same way as the one which we discussed above. Rather 
the correct representation corresponds to an alternative embedding in which: 

50(8)i ^ 50(7)i X Is 

8 ^{8, a). (6.4) 

Here, the 8 on the right-hand-side of this equation is the 8-dimensional spinor representation 
of 50(7). [ See any book on Lie groups, e.g. Ref.( ^8]), regarding the properties of the 
orthogonal groups.] In general, the scaling dimension, x, of a Kac-Moody primary field 
transforming under a representation p of a group is given bylll: 

- - (6.5) 



CA + k' 



where Cp is the quadratic Casimir invariant for the representation p, Ca is the quadratic 
Casimir for the adjoint representation and k is the Kac-Moody central extension (level). For 
the spinor (8), vector (7) and adjoint (21) representations of 5*0(7), the Casimirs are: 

O7 = 3 
Cs = 21/8 

0*21 = 5. (6.6) 
Taking A; = 1, we obtain for the scaling dimensions 

xr = 1/2 

X8 = 7/16. (6.7) 
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We see that a consistent conformal embedding represents the 8 Majorana fermions, xa, 
of S0{8)i, in terms of the S0{7)i spinor field, 5"^, {A = 1,2, 3, ...8.) and the Ising order 
parameter field, a: 

Xaoc Sa- (J- (6.8) 

Since = 1/16, the dimension adds up correctly to 1/2 for a free fermion. Comparing to 
our previous bosonization formula of Eq. ( |3.28| ) we note that the Ising factor, a, is common 
to both formulae and that the SO (7) spinor field has replaced the product of SU{2) spin 
and isospin fields. The free fermion spectrum of Table I can be rewritten in terms of 
5*0(7)1 X Ising conformal towers. We may now use the conformal embedding of 5*0(7)1 
conformal towers [labeled by dimensions] into 5f/*^(2)i x 5f/*^(2)i xSU^ {2)2 conformal 
towers [labeled by angular momenta]: 

7^(l/2,l/2;0) + (0,0;l) 

8 ^(1/2,0; 1/2) + (0,1/2; 1/2), (6.9) 

(Note that both KM algebras have the same conformal anomaly parameter: c = 7/2 = 
1 + 1 + 3/2.) From the equation above and TABLE I, we find for the combinations of 
conformal towers of the free fermion spectrum: (50(7)i, Ising) = (1, 1), (8, a), (7, e). 

It is now immediately evident that the 5*0(7) symmetry will be preserved in any new 
spectrum obtained from this one by fusion with an Ising field only. In particular, the 
combinations of conformal towers at the non-trivial fixed point are given in this notation 
by: (5*0(7)1, Ising) = {^,<^),{8,l),{8,e),{7,a). Note that the two-fold degenerate sets of 
representations in Table III at x — 1/16 = 3/8, 1/2 and 7/8 collapse to single representations 
of 5*0(7) in all three cases. The conformal tower at x — 1/16 = 1 in Table III actually 
corresponds to a descendent of (1,0") from the 50(7) point of view. We may also simplify 
Table V by constructing the spectrum up to s — 1/16 < 2 in this new basis. The first 
descendents of the 50(7) conformal towers are given by: 

5*0(7)1— primary — ^ 5*0(7)i —descendant 
1 ^ 21' 

7 ^ 7', 35' 

8 ^8', 48'. (6.10) 

We find that all states in Table V can be grouped into 5*0(7) representations and this 
explains all "accidental" degeneracies. (Of course other degeneracies persist which arise 
from the fact that the spacing of levels of all conformal towers in any conformal field theory 
is Sx = 1.) 

We remark that there is an isomorphism^ of the fusion rules of the 50(7)i and the 
Ising conformal field theories, corresponding to the following mapping of 5*0 (7) 1 and Ising 
conformal towers: 

5*0(7)1 Ising 
1 ^ 1 

8^ a 

7^e (6.11) 
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Note that since the free fermion spectrum corresponds to a sum of products of the three 
pairs of identified conformal towers, namely {S0{7)i, Is) = (1,1), (8, a), (7, e), we can go 
from the free fermion to non-trivial fixed point by fusion either with a in the Ising sector 
or, equivalently, with 8 in the SO {7)i sector. 

The decomposition of the adjoint representation of 5*0(7), under which the 5*0(7)1 
current [a KM descendant of the 5*0 (7) i identity operator] transforms, gives information 
about the additional conserved quantities: 

50(7)i ^ SU''{2)i X SU'^{2)i x SW{2)2 

21' ^ (1', 0; 0) + (0, 1'; 0) + (0, 0; 1') + (1/2, 1/2; 1). (6.12) 

At the zero coupling, infinite coupling and non-trivial fixed points, there are 21 (Hermitean) 
conserved quantities with these quantum numbers. The first 9 clearly correspond to the 
currents, Ji, I2 and J. The remaining 12 operators are uniquely identified by their quantum 
numbers as ipiailj2, ipi(j'^Sip2 together with the hermitean conjugates of these operators. 
We see that one linear combination of these corresponds to the extra conserved quantity 
identified previously in the NRG work, Eq. (|6.1|) . 

The 28 fermion bilinears that can be formed from 8 Majorana fermions, which transform 
under the adjoint representation of 50(8) and represent the 5*0 (8)1 KM current operator, 
decompose into the 21 and 7 dimensional representation of 50(7). Thus under the conformal 
embedding 5*0(8)i 5*0(7)i x Js, 21 of the fermion bilinears become the 5*0(7) current and 
the remaining 7 become the product of the vector (7) and e primaries. It can be seen that 
both types of bilinears appear in the Kondo interactions. As stated above, these interactions 
do not manifestly preserve the 5*0(7) symmmetry for any non-zero value of the couplings; 
rather this symmetry is "dynamically restored" at the non-trivial critical point. 

This (non-abelian) bosonization scheme of the 4 species of Dirac Fermions in terms of 



an 50 (7)1 Wess-Zumino-Witten theory and an Ising model [first pointed out in Ref. ( |Tj) ], 
appears to be closely related to a physical picture for the 2-impurity Kondo problem proposed 
more recently by Sire et al.El and Gan.0 In these papers, the fermions are first transformed 
by ordinary abelian bosonization. Certain linear combinations of bosons are defined and 
then it is observed that a subset of the boundary operators can be re-expressed in terms of 
some different fermions, a transformation sometimes referred to as "refermionization" . For a 
special, anisotropic, value of the Kondo couplings, up to irrelevant operators, it is found that 
only 1 of the new Majorana fermion couples to the impurity spins, the other 7 completely 
decoupling. This theory has an obvious 5*0(7) symmetry. This procedure apparently gives 
an explicit realization of the 5*0(7) bosonization scheme, with the added bonus that the 
5*0(7) symmetry at the non-trivial critical point becomes manifest. 



VII. DISCREPANCY WITH OTHER CALCULATIONS 

Finite temperature Monte Carlo (MC) calculations by Fye and Hirsch.0 have so far 
discovered no evidence for the non-trivial critical point discussed here. It was suggested 
that the non-trivial critical point may be an artifact of the energy-independent coupling 
constant approximation used in the NRG work. Furthermore, NRG calculations by Sakai et 
al.0 on the related 2-impurity Anderson model, while reproducing the results of Ref. ( ^ in 
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the energy (and parity) independent coupling constant case, did not find such a critical point 
using energy-dependent coupling constants or a "parity-splitting" term. In this section we 
discuss the apparent discrepancy between these calculations and the ones presented in the 
present paper. 

Our most important point was derived in detail in Sec. II and Sec. V; it is the existence of 
two different types of P-H symmetry, only one of which leads to a non-trivial fixed point. As 
explained in Sec. II, when the first type of particle-hole symmetry is maintained, some sort of 
phase transition as a function of inter-impurity coupling, K, is inevitable.! The reason is that 
for large antiferromagnetic K there is no phase shift and for large ferromagnetic K the phase 
shift is tt/2 in both channels. Furthermore, because of this type of particle-hole symmetry, 
the phase shifts, if they are well-defined, can only be or tt/2. Both of these phases are 
absolutely stable (no relevant or marginal operators). Thus they must be separated by (at 
least) one point which is not in either phase. In principle, this could either be a first order 
transition point or a critical point. This argument seems to us very general and convincing. 
It should not, for instance, depend on ignoring energy- dependence of the coupling constants, 
discussed in Sec. II. It only relies on the rather well-established local Fermi liquid picture of 
the single-impurity Kondo effect developed by Nozieres and others. The statement that the 
phase shift is 7r/2 for some range of parameters has a precise mathematical meaning. It is 
a statement about the asymptotic behaviour of the electron Green's function, discussed in 
Sec. V. If the Green's function did not have this behaviour for a range of ferromagnetic K 
this would be rather shocking and would presumably indicate a breakdown of the standard 
theory of the single-impurity Kondo effect. Likewise for a range of antiferromagnetic K. We 
think that measuring the quantity, S'(i) , governing the asymptotic behaviour of the electron 
Green's function, defined in the previous section, would be the easiest way to verify the 
existence of the non-trivial critical point. It ought to be possible to show that S'(i) = 1 
for sufficiently large antiferromagetic inter-impurity coupling and 5(1) = — 1 for sufficiently 
large ferromagnetic inter-impurity coupling. A numerical study of how and where it passes 
between 1 and —1 would give information about the non-trivial critical point. We emphasize 
that 5(1) 7^ ±1 corresponds to non- Fermi liquid behaviour. At T = some sort of phase 
transition must exist between these two Fermi liquid cases. An analysis of scaling with 
temperature will be required to study this transition. 

Of course, this argument says nothing about the nature of the phase transition separating 
the two stable phases. It could, for example, be first order, in which case there would not 
be a critical point. To determine the nature of the phase transition requires other methods. 
The NRG and CFT methods provide complementary approaches, which taken together, 
suggest rather clearly that the two phases are separated by the non-trivial critical point 
discussed in this paper. The CFT approach allows a systematic classification of possible 
critical points, given the assumption of conformal invariance at the critical point and the 
various symmetries of the problem. The NRG results allow a determination of which of 
the possible critical points actually occurs. Having identified a candidate critical point, and 
shown that it does occur in a particular microscopic formulation of the model (corresponding 
to the NRG with energy-independent coupling constants), we can then determine whether 
it will be stable under changes in the microscopic Hamiltonian using the RG. The complete 
set of relevant and marginal operators at the proposed critical point are those in Table VI, 
together with the marginal current operators, J, If the microscopic Hamiltonian has the 
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first type of particle-liole symmetry (togetlier witli parity and at least a U{1) subgroup of the 
spin-rotation symmetry) then all but one of these operators are forbidden to appear in the 
effective Hamiltonian of Eq. ( p.4| ). The presence of the first type of particle- hole symmetry 
impliesiil@ the diagonal isospin symmetry, /. This diagonal isospin symmetry together 
with spin symmetry [or at least a U{1) subgroup of it] forbids all these operators except e 
and an element of the eighth operator in Table VI, which is forbidden by parity. (See the 
beginning of Sec. V.) Thus we reach the crucial conclusion that the non-trivial critical point 
can be reached by adjusting only one parameter, the coupling constant corresponding to e in 
the effective Hamiltonian, assuming the first type of particle-hole symmetry. By varying the 
inter-impurity coupling, K, or any other parameter, so as to pass from the inter-impurity 
singlet phase to the Kondo-screened phase, we can make this coupling constant pass through 
zero. Thus this is a conventional critical point; it is reached by varying only one parameter, 
assuming the first type of particle-hole symmetry. 

This analysis indicates that allowing energy-dependent coupling constants should not 
make any difference, provided that the first type of particle-hole symmetry is maintained. 
No lowering of symmetry ensues from allowing energy-dependent coupling constants, as long 
as the first type of P-H symmetry is maintained, so no additional operators are permitted 
in the effective Hamiltonian describing the critical point. Of course, the actual location of 
the critical point, as a function of the various Kondo and RKKY couplings can change. 

The NRG work on the Anderson mo del0 considers several different Hamiltonians. The 
model with energy- independent coupling constants, reduces, in the Kondo limit of the An- 
derson model, to the same model studied previously^; in this case apparently the same 
non-trivial critical point is observed. The other models considered, with energy-dependent 
coupling constants or a "parity-splitting" term do not have the first type of particle-hole 
symmetry. Thus, by the discussion, in Sec. II, we do not expect them to exhibit the non- 
trivial critical point, and indeed they do not. Some of these models do have the second 
type of particle-hole symmetry. This symmetry is not enough preserve the critical point as 
we discussed in Sec. II, V. However, it should be possible to study models with energy- 
dependent coupling constants which do not break the first type of particle-hole symmetry, 
and which therefore should exhibit the non-trivial fixed point. This corresponds to choosing 
couplings Ne^o{E) [corresponding essentially to Wq^i{E) in the notation of Ref. ( |l3)] which 
are non-trivial even functions of energy, E, for both parities, [p = 0, 1]. (The other param- 
eters in the Hamiltonian must also be chosen to preserve particle-hole symmetry.) Further 
calculations with such models would be valuable in settling this controversy. 

Most of the QMC work of Ref. ( 0) was done on the one- dimensional tight-binding model 
at half-filling with 5-function Kondo interactions, of the form: 

J['ilj\R/2)aij{R/2) ■ Si + ^\-R/2)a^{-R/2) ■ S^s], (7.1) 

where R = na, a being the lattice spacing. [In the case n odd, the origin is midway between 
2 sites.] The various cases, = 1,2,4, 8 were studied. We see from Eq. ( p.S] ) that the first 
type of P-H symmetry occurs for n even and the second for n odd. A systematic search for 
the non-trivial critical point was only made in the case n = 1, where it is not expected to 
occur! The reason for this was connected with the fact that these authors did not put in an 
inter-impurity coupling by hand. They expect (based on weak-coupling perturbation theory) 
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that the RKKY couphng will be ferromagnetic for n even and antiferromagnetic for n odd. 
The NRG work of Ref. ( suggests that the non-trivial critical point should occur for 
antiferromagnetic inter- impurity coupling, of O(T^). From another point of view, it is quite 
easy to see that these models are generally in the Kondo- screening phase at strong coupling, 
so if a transition is to be encountered they had better be in the inter-impurity singlet phase 
at weak coupling. A necessary condition for this is presumably an antiferromagnetic RKKY 
coupling. 

There is another potential problem with the even n models. In this case Nf.{0) = 0, 
[see Eqs. ( ^.11 ), so the even channel has only irrelevant Kondo couplings. If we 



assume the even channel decouples, we obtain a single- channel Kondo effect with a spin-1 
impurity for ferromagnetic inter-impurity coupling. This leads to a Fermi liquid fixed point, 
with an underscreened s = 1/2 effective impurity. However, we expect that higher orders 
of perturbation theory will generate a non-zero Kondo couplings to the even channel at 
E = 0. Such a coupling to the leftover s = 1/2 impurity would be relevant, provided that 
it is antiferromagnetic. However, if it is ferromagnetic, the underscreened Fermi liquid fixed 
point would be stable and the symmetry argument given above for the existence of the 
non-trivial critical point may fail. 

Thus it seems quite likely that the non-trivial critical point will not occur for any value of 
n and J in this model. However, it should be possible to see the critical point by generalizing 
the model somewhat, while maintaining the desired symmetry. It may be sufficient to take 
an even n model (say n = 2) and add a direct inter-impurity coupling, K. To avoid the 
potential difficulty mentioned in the previous paragraph, it may also be necessary to modify 
the Kondo couplings somewhat while preserving the first type of P-H symmetry. It can be 
seen that this symmetry occurs at half-filling provided that there is site-parity, ie. refiection 
symmetry about a site as occurs for n even. The second type occurs if there is link-parity. 



ie. reflection about a link. One example would be to make the replacement in Eq. (jTA 



^{-R/2) ^Votlj{0) + iVitlj{-a). (7.2) 
This gives, in wave-vector space: 

v{k) = vo + ivie'^'K (7.3) 

The first type of P-H symmetry holds provided that the fj's are real. For vq and Vi both 
non-zero, A'^e(O) and No{E) are both non-zero. However, we find that the RKKY coupling 
is always ferromagnetic for all the models of this type that we have considered with the first 
type of P-H symmetry. Thus it may be necessary to add a direct inter-impurity coupling 
which is varied to pass between the two stable phases. 

It may still be difficult to see the non-trivial critical point by studying the staggered 
susceptibility since it only shows a logarithmic divergence. It may be necessary to choose 
couplings very close to the critical value and very low temperatures to see the anomalous 
behavior. As discussed above, the single-particle Green's function would likely exhibit a 
much clearer signal of the critical point. 

In conclusion, we expect that further numerical work, along the lines of Refs. ( |10|) and 
( 1^) should be able to see the critical point provided that: 
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1) the first type of particle-hole symmetry is maintained, 

2) a parameter is varied to pass from the inter-impurity singlet to Kondo screened phase 
and 

3) sufficiently low temperatures and finely tuned parameters are obtained. 

The last condition should be much easier to achieve if the single-particle Green's function is 
measured instead of the staggered susceptibility. 
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TABLES 



TABLE L 
(3^). 


Conformal towers 


occurring for free fermions with the boundary conditions of Eq. 


h 


i2 


j 


Ising 


X 




1/2 





1/2 
1/2 





1/2 


1/2 
1/2 




1/2 

1/2 

1 

1 




1 

a 
a 

€ 

1 

e 




1/2 

1/2 

1 

1 

1 


TABLE IL Conformal towers 
to 4) [or e or (1/2,1/2)] fusion. 


occurring for free fermions with a 7r/2 phase 


shift, corresponding 


h 




j 


Ising 


X 




1/2 





1/2 
1/2 





1/2 


1/2 
1/2 


1 

1/2 
1/2 



1 


1 
a 
a 
e 
1 
e 


1/2 
1/2 
1/2 
1/2 
1/2 
3/2 




TABLE IIL 


Conformal towers occurring 


after a fusion. 




h 




J 


Ising 


x 1 
16 




1/2 





1/2 
1/2 


1/2 





1/2 


1/2 


1/2 
1/2 




1/2 
1/2 
1 


1/2 
1/2 
1 


a 
1 
1 

a 
a 

e 
e 
a 




3/8 
3/8 
1/2 
1/2 
7/8 
7/8 
1 
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TABLE IV. All first descendents of all conformal towers. The subscript in the first row labels 
the Kac- Moody level, k. SU{2) descendents are labelled by their spin. All descendents will be 
marked '. 



(0)i 


a), 


(0)2 




(1)2 


1 


a 


e 


V 


2 


1' 


1' 3' 
2 ■ 2 


0'. 1' 




a' 
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TABLE V. Spectrum at the non-trivial fixed point. 





j 


Ising 


X 16 




%RG 


(0,0) 





(T 





0+ 


.00000 


(io), (o,i) 


1 

2 


1 


3 
8 


1+ 
2 

I — 
2 


.37761 
.37764 


(0,0) 


1 


a 


I 
2 


0~ 


.50454 


V2' 2/ 





a 


1 
2 


1+ 


.50724 
.50726 


(io), (o,i) 


1 

2 


e 


7 
8 


1- 
2 

1 

2 


.88696 

.88702 


(0,0) 


1' 


a 


1 


0+ 


.99952 


V2' 2/ 


1 


a 


1 


0+ 
1~ 


1.00296 
1.00298 


(l',0),(0,l') 





a 


1 


1- 

1 + 


1.00642 

1.UUD44 


(0,0) 





a' 


1 


0+ 


1.01078 


(io), (o,i) 


2 


1 




1+ 

2 

1- 
2 


1.39056 
1.39066 


(2 'U),(U, 2 j 


1 

2 


1 


1 3 

Is 


1 + 

2 

1- 
2 


i.oyozi 

1.39329 


(1 > 2i'i2'-'- ) 


1 
2 


1 




3- 
2 

1 + 
2 

1- 
2 

3 + 
2 


1 QQRQQ 

i.oyoyy 
1.39700 
1.39709 
1.39716 


(i0),(0,i) 


2 


1 


-^8 


1 + 
2 

1 ~ 

2 


1.52680 
1.52687 


(0,0) 


1' 


(7 


ll 




1.50437 


(1' 1) (1 I') 





(T 




0" 

1+ 
1- 

0+ 


1.51090 
1.51094 
1.57554 
1.57556 


(0,0) 


0' 


a 


1^ 


Q- 


1.56562 


V2' 2/ 


1' 


a 




0- 

1+ 


1.56892 
1.56895 


(l',0),(0,l') 


1 


a 




1+ 
1~ 


1.57224 
1.57226 


(0,0) 


1 

i 


(T 


1 1 


U 


1 tiT A ceo 


V2' 2/ 





a' 




Q- 

1+ 


1.62305 
1.62312 


(io),(o,i) 


1 

2 


e' 


i| 


1- 
2 

1 + 
2 


1.92017 
1.92029 


(o,|),(io) 


3' 
2 


e 




1- 
2 


1.97512 
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1+ 
2 


1.97514 


(2 2 ) 


1 

2 


1 7 

e Ig 


1- 
2 


i.y /oi / 








1+ 
2 


1 Q7fi9/1 


/I -\/\ r-i/ l\ 
(2'1 '2) 


1 

2 


e Ig 


1- 
2 


i.yozz4 






3+ 
2 


1.98225 








1+ 
2 


1.98230 








3- 
2 


1.98236 


/I rl^ /n 1 \ 


1' 

2 


1| 


1- 
2 


z. 23938 








1+ 
2 


2 23Q47 




TABLE VI. Operator content at the non- 


-liiviai ciiLieai poiiiL. 




«1 




i 


Ising 













1 





1/2 





1/2 


a 


1/2 


1/2 





1/2 


a 


1/2 





1/2 


1/2 


a 


1/2 





1/2 


1/2 


a 


1/2 








1 


1 


1/2 











e 


1/2 


1 /9 


1/2 







1 /9 








1 


e 


1 


1/2 


1/2 


1 


1 


1 


1/2 


1/2 





e 


1 


1/2 


1/2 


1 


e 


3/2 
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